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Portugal’s leadership in sustainable energy and to attain the ambitious goals set in the government 
programme, namely consolidating the renewable energies cluster in Portugal, which will represent 
approximately by 2020 more than three times the 35,000 jobs estimated in 2010 and further 
developing the industrial cluster related with energy efficiency, creating 21,000 new jobs. 

Keywords: The main purpose of this paper is to perform a prospective study and to discuss the various factors 
Employment effects that influence the appraisal of sustainable systems integrating environmental, social and economic 


Input-Output analysis . i dimensions, mainly focusing on the RES-E jobs case by means of a multiobjective Input-Output model, 
Multiobjective interval linear programming 


(MOLP) based on current data availability. 
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1. Introduction 


The benefits of renewable energy are widely accepted by 

fostering a decrease in Green House Gases (GHG) emissions and 

ee an improvement in energy security. Moreover, it is commonly felt 

* Corresponding author: INESC Coimbra, Rua Antero de Quental 199, 3030-030 that the move towards a greener economy will create a large 
Coimbra, Portugal. Tel.: +351 239 851040; fax: +351 239 824692. l & y : jargi 

E-mail addresses: coliv@inescc.pt (C. Oliveira), dcoelho@isec.pt (D. Coelho), number of new jobs across many sectors and that this fact will act 

patsilva@fe.uc.pt (P. Pereira da Silva), cantunes@inescc.pt (C.H. Antunes). as a vital incentive to sustainable development and to economic 


1364-0321/$-see front matter © 2013 Elsevier Ltd. All rights reserved. 
http://dx.doi.org/10.1016/j.rser.2013.01.011 


C. Oliveira et al. / Renewable and Sustainable Energy Reviews 21 (2013) 444-455 445 


growth. The exploitation of renewable energy sources in electri- 
city generation is expected to have a crucial contribution to the 
overall rise of employment opportunities in several sectors, 
including equipment manufacturing, construction, administrative 
and service activities [29], but this claim stands on an uncertain 
footing. While studies often present renewable energy as a boost 
to the economy through the generation of large job growth (such 
as [3]) there are often overly optimistic or simplistic assumptions 
that lead to these results. In some cases reasonable assumptions 
are made, but selective reporting of results can lead to a false 
impression of job creation. 

Studies analyzing the impact of job creation in the RES-E sector 
can be broadly categorized into two main categories: input- 
output (I-O) methods and analytical methods, both with their 
own distinct advantages and disadvantages [16, 30]. While I-O 
methods are easily used within a National scope, analytical 
methods are commonly used for regional or provincial studies 
(e.g., [3, 27]. On the other hand, since analytical methods usually 
account for direct employment effects only, traditional methods 
applied in jobs and economic impacts analysis mostly rely 
on I-O models to estimate employment creation or loss. This kind 
of models allows for the representation of the economy as a 
whole, recording the flows of goods and services industries 
trade with one another. Those flows are registered in an I-O 
matrix, simultaneously by origin and by destination which illus- 
trates the relationship between producers and consumers as well 
as interdependencies of industries for a given year [17]. Therefore, 
these models allow capturing employment multiplier effects 
and macroeconomic impacts of shifts between sectors, account- 
ing for losses in one sector created by the growth of another 
sector [5]. 

Many countries compile I-O tables for their economies at 
regular time intervals both as a national statistical requirement 
and for the purpose of providing detailed databases for policy 
analysis. Therefore, I-O analysis is an analytical tool adequate for 
the evaluation of the inter-relations between different economic 
activities being often applied to assess economy-energy-envir- 
onment (E3) interactions [12]. 

I-O analysis limitations mainly refer to the hypotheses 
assumed within the model. In fact, this modelling technique is 
based on a set of assumptions that might be considered as its 
main drawbacks. The technical coefficients are considered as 
constant over time, there are no economies or diseconomies of 
scale in production or factor substitution and it is also assumed 
homogeneous output for each activity sector. Another issue is 
that, in general, final demand is exogenously determined. Finally, 
there is the assumption that there are no bounds for the capacity 
of production, that is, supply is supposedly infinite and perfectly 
elastic. 

Albeit apparently these restraining assumptions limit the 
application of I-O analysis, it is possible to eliminate or avoid 
some of them through adequate adjustments. For instance, the 
uncertainty handling of the technical coefficients may be tackled 
through the use of interval or fuzzy programming techniques 
(see [23] and [4], respectively). 

The use of I-O multipliers is particularly appropriate for 
evaluating the contribution of a particular industry (e.g., RES-E) 
to the economy and for performing the impact assessment of 
broad policy instruments. However, in order to get the whole 
picture of the impacts of an economic activity and an environ- 
mental impact assessment, a multi-objective analysis must be 
implemented as well [6]. Therefore, an overall analysis of the 
trade-offs regarding the conflicting axes of evaluation intrinsic to 
sustainable development will be performed by means of a multi- 
objective I-O model with interval coefficients, based on current 
data availability. 


In the next section of this paper, the methodological frame- 
work used herein will be briefly presented, followed by the 
explanation of the implementation of the methodology in Portu- 
gal. Then, some illustrative results will be analyzed, considering 
different sources of uncertainty and scenarios and, finally, some 
final remarks and future work developments are drawn. 


2. I-O model description 


I-O matrices allow the representation of each sector’s produc- 
tion process through a vector of structural coefficients that 
describes the relationship between the intermediate inputs con- 
sumed in the production process and the total output. The supply 
side is split into several processing industries that deliver their 
total output (production), for intermediate consumption or final 
demand. These relationships can be illustrated through the 
following equation: 


Xj = So xj+yi (2.1) 


where x; is the output of sector i, xj is the input from sector i to 
sector j, and y; is the total final demand for sector i. 

The monetary values in the transactions matrices can then be 
converted into ratios called technical coefficients. This is done by 
dividing each cell of the domestic intermediate matrix by its 
column total (output at basic prices). 

Considering the hypothesis of constant returns to scale, 
Eq. (2.1) becomes: 


n 
x= 5 aijXj+Yi (2.2) 
f=1 


in which the coefficients aj are the amount of input delivered by 
sector i to sector j per unit of sector’s j output, known as 
technological coefficients (or direct coefficients). 

The productive system at a national level can then be repre- 
sented through the following basic I-O system of equations: 


X=Ax+y (2.3) 


where A is a matrix of technological coefficients, y is a vector of 
final demand, and x is a vector of the corresponding outputs. 

In order to finally calculate the output multipliers, one needs 
to derive Leontief inverse matrices. 

Eq. (2.3) can then be rearranged to 


x= (I-A)"'y, (2.4) 


where I is the identity matrix with convenient dimensions and 
(I—A)~' is also known as the Leontief inverse. Each generic 
element, bj, of (I—A)~’ represents the total amount directly and 
indirectly needed of good or service i to deliver one unit of final 
demand of good or service j. 

Several empirical contributions have used the seminal meth- 
odological approach developed by Leontief. For instance, the 
evaluation of the inter-relations between different economic 
activities established by I-O analysis may be used to assess 
economy-energy-—environment interactions—see, for example, 
the work developed by Gay and Proops [11], Peet [25] and Mu 
et al. [19]. Usually, these empirical applications of I-O analysis 
highlight the use of the Leontief inverse indicating the direct and 
indirect requirements of production that are needed to satisfy a 
particular final demand vector, being also known as the multiplier 
matrix. 
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3. Employment multiplier concepts 


Although precise definitions vary, direct jobs are related 
to a sector’s core activities, such as feedstock conversion, manu- 
facturing, project development (including site preparation and 
installation) and operations and maintenance of the different 
components of the technology, or power plant, under considera- 
tion [30]. 

Therefore, the direct contribution of an industry in terms of 
output or employment can easily be measured by its level of 
output or the number of workers in the sector, respectively. Since 
the employment to output ratio is given for each sector in an I-O 
table, the overall significance and contribution of an industry to 
total employment can also be calculated by assuming that the 
sectorial employment ratios are fixed. 

Indirect jobs refer to the “supplier effect” of upstream and 
downstream suppliers, corresponding to the industrial input 
sectors in the production and the operation and maintenance of 
renewable energy technologies. Examples might include the jobs 
required to extract and process raw materials, such as steel for 
wind turbine towers as well as positions in government minis- 
tries, regulatory bodies, consultancy firms and research organiza- 
tions working on renewables [30]. 

Thus, the indirect contribution of an industry to either total 
output or employment is not simply observable unless the 
multiplier and flow-on effects are taken into account. Therefore, 
the employment multiplier may be interpreted as the impact on 
the overall employment if the final demand in sector j increases 
by one unit. The employment multiplier for sector j, Eʻ, is thus 
defined as follows: 


n 
EP = 5 eibij, 


i=1 


(3.1) 


where e; denotes the number of persons with full time employ- 
ment per one Euro output for each sector i, by is the (i, j)th 
element of the closed Leontief inverse matrix and n is the number 
of sectors. These multipliers would represent the number of new 
jobs created expressed as total employment for every new 
employee to meet increased final demand of new output, but 
one may wish to relate the simple or total employment effect to 
an initial change in employment, not to final demand (and 
output) in monetary terms. In this situation the employment 
multiplier, Ej, is: 


(3.2) 


4. The MOLP I-O model with interval coefficients 


I-O analysis and linear programming (LP) are closely related. 
In its simplest form, I-O analysis may be regarded as a simple 
particular case of LP [8]. The use of the I-O methodology in the 
framework of LP models allows obtaining value-added informa- 
tion, which would not be possible to achieve with the separate 
use of both techniques. Inter/intra-sector relations embedded in 
I-O analysis allow designing the production possibility frontier. 
LP models enable choosing the optimum level of activities to 
optimize a given objective function, satisfying the production 
sector relations imposed by I-O analysis. Traditional studies, 
which use I-O analysis in the framework of LP, generally consider 
a single objective function, usually an aggregate economic indi- 
cator. However, in most real-world problems multiple, conflicting 
and incommensurate axes of evaluation of the merit of potential 
solutions are inherently at stake. In this context, mathematical 


programming models for decision support become more repre- 
sentative of reality if distinct aspects of evaluation are explicitly 
considered. In this context, an MOLP I-O model to deal with E3 
interactions has been proposed elsewhere [23, 24]. Some changes 
are now incorporated into the model herein proposed: an 
updated data set of I-O symmetrical product by product tables 
for the total flows at basic and current prices of 2008, the 
construction of distinct vectors for the incorporation of RES-E 
production within the I-O matrix, the use of multipliers to obtain 
the impacts of RES-E on the overall employment and the con- 
sideration of different objective functions more consistent with 
the aim of this study (the maximization of Gross Domestic 
Product—GDP as a proxy for economic growth, the maximization 
of the overall employment in the economy and the maximization 
of RES-E production). 

The model includes two main types of constraints: coherence 
constraints (based on I-O analysis) and defining constraints. The 
economic and environmental defining constraints with interval 
coefficients have been imposed with interval (upper/lower) 
bounds consistent with data available [14,18]. 


4.1. Model constraints 


The matrices are given in capital letters, the vectors in small 
bold letters, the scalar elements are given as small letters and the 
letter T designates the transpose of a matrix or a vector. 


4.1.1. Coherence constraints 

The intermediate consumption and final demand of goods or 
services of each activity sector shall not exceed the total amount 
available from national production and competitive imports of 
that same good or service. 

Energy and non-energy branches: 


Ax+ Dy <imp‘+x (4.1) 


A is the technical coefficients (product-by-product) matrix, where 
each element is the amount of good or service i needed to produce a 
unit of good or service j. This matrix has interval coefficients and 
is given in hybrid units (toe/million of Euros or million of Euros/toe). 
D is the diagonal matrix whose main diagonal corresponds to the 
coefficients of consumption of energy and non-energy commodities 
by the final demand sectors (households and non-profit institutions 
serving households—NPISH, public consumption, gross fixed capital 
formation, acquisitions less disposals of valuables and stock changes 
and exports); 


x is the national output vector. 

y is the final demand vector. 

imp‘is the vector of competitive imports with no endogenous 
substitutes (for energy commodities only). 


I-O energy consumption coefficients are given as intervals to 
cope with their uncertainty. The upper and lower bounds of these 
intervals correspond, respectively, to a more or less pessimistic 
stance, regarding the energy coefficient settings for the planning 
horizon (2020). The coefficient set leading to the broadest feasible 
region incorporates an improvement of the energy efficiency 
measured through the energy use coefficients of the different 
activity branches. It takes into consideration a reduction of 5% for 
the electricity, diesel oil and gasoline consumption coefficients. 
The energy coefficient setting leading to the more stringent 
version of the feasible region considers an increase of 5% in 
natural gas consumption coefficients regarding co-generation 
and thermoelectricity generation, since national authorities are 
replacing fuel oil consumption by natural gas. 
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4.1.2. Defining constraints 


4.1.2.1. Economic constraints. Several consumption representations 
are considered in the model: the households’ consumption in the 
territory (consumption in the territory by resident and non-resident 
households), the residents’ (households and NPISH) consumption, 
the resident households’ domestic consumption which is linearly 
dependent on the available income, and the tourism imports given 
as a proportion of the resident household consumption [21]. 

With respect to foreign trade, it is possible to obtain: total 
exports (excluding and including tourism) at constant FOB (free 
on board) prices, exports at constant purchasers’ prices, total 
exports (excluding tourism) at CIF (cost, insurance and freight) 
prices, total imports (excluding tourism) at CIF prices, and total 
imports (excluding and including tourism) at FOB constant prices. 

Gross domestic product (GDP) is computed according to the 
production approach (gdpproq) and the expenditure approach 
(gdp). 

GDP production approach: 


EAP prod = Zav + ts (4.2) 


where gav is the gross added value (see (4.3)) and ts is the sum of 
the components of vector of net taxes ts (see (4.4)). 
Gross added value (gav): 


gav =a! x+a!x—alx+al.x (4.3) 


where aw, at, as, aos are the vectors with the proportion of wages, 
taxes, subsidies, gross operating surplus and gross mixed incomes 
on the total output of each branch. 

Net taxes (ts): 


ts = AtsX + Dtsy (4.4) 


where Ats is the matrix with the proportion of net taxes on goods 
and services on the total output of each branch, and Des is the 
diagonal matrix whose main diagonal is the vector with the 
percentages of net taxes on goods and services aimed at house- 
hold/NPISH/public consumption/gross fixed capital formation/ 
changes in inventories/acquisitions less disposals of valuables/ 
exports on their respective total values. 
GDP expenditure approach 


gdp =rc+g+efcf+sc+aldv+expfob—mfob (4.5) 


where rc is the resident’s consumption, g is the public consump- 
tion, gfcf is the gross fixed capital formation, sc is changes in 
inventories, aldv is acquisitions less disposals of valuables, expfob 
is exports at FOB purchasers’ prices (including tourism) and mfob 
is imports at FOB prices (including tourism). 

The GDP at current prices is obtained from the distinct 
components of GDP (expenditure approach) at constant prices, 
which are multiplied by the corresponding deflators. 

The residents’ disposable income is equal to the difference 
between the National Available Income and the sum of the 
available income of corporations and public administration. 

The employment level (emp) is obtained by using labour gross 
productivity coefficients (17) for each branch: 


emp = Ix (4.6) 


Energy imports (energy external dependence) are obtained in 
the following way: 


(e:)"imp* + (e2) (Ap x+ Dirty) (4.7) 


where e; and e, are vectors of ones with convenient dimensions, 
An is the matrix of non-competitive import coefficients of energy 
(given as intervals to capture uncertainty) and Di" is the diagonal 
matrix whose main diagonal is the vector with the percentages of 


non-competitive imports of energy goods aimed at household/ 
NPISH/public consumption/gross fixed capital formation/changes 
in inventories/acquisitions less disposals of valuables/exports on 
their respective total values (also given as interval coefficients). 


4.1.2.2. Environmental constraints. CO, emissions from fuel 
combustion are easily obtained from the I-O table, where the 
total fuel use is the total amount of fuel production plus imports. 
Nevertheless, the energy use for exports and investment shall not 
be taken into account in the emission computations [26]. 


coz = DyjD-DolAcx + Dey—Ney Nex](4$) (10°) (4.8) 
where Diy is the diagonal matrix with conversion factors from toe 
to terajoules (TJ) for each type of energy, Dcis the diagonal matrix, 
whose main elements are the carbon (C) emission factors for each 
type of energy (given as interval coefficients), Dois the diagonal 
matrix, whose main elements are the fractions of carbon oxidized 
for each type of energy, Ag is a sub-matrix of matrix A with the 
energy consumption coefficients, Dg is a sub-matrix of matrix D 
with the energy consumption coefficients, Ng is the matrix with 
the coefficients of energy use with non-energy purpose for each 
branch and 44/12 is the ratio between the molecular weights of 
CO, and C. 


4.2. Objective functions 


The allocation of energy resources shall be made having in 
mind that the energy sector is a part of the economic system as a 
whole and that energy planning requires the consideration of 
economic, social, energy and environmental objectives, the four 
main pillars of sustainable development. In this way, the model 
herein proposed considers the following objective functions: (a) 
Maximization of GDP (economic growth); (b) Maximization of the 
level of employment (social welfare); (c) Maximization of RES-E 
production (energy and environmental concerns). 


5. Implementation of the methodology in Portugal 


The I-O symmetrical product by product tables for total flows 
at basic prices at current prices of 2008 herein used were 
produced by the Portuguese Department of Foresight and Plan- 
ning and consist of 64 production sectors [7]. A hybrid I-O unit’s 
matrix has also been adjusted for the base year of this study based 
on the Portuguese Energy Balance and on the energy statistics 
available at the Portuguese Directorate General of Energy and 
Geology for the base year of this study (2008), although a 
monetary unit’s matrix has been used to compute the RES-E 
employment multipliers. 

The most straightforward way to assess employment effects is 
to collect employment data directly from the units involved in the 
considered activities. However, in most countries, such as Portu- 
gal, the number of other related RES-E activities are still few and 
reliable data regarding employment inputs are either scarce, 
missing or vary significantly from case to case. Therefore, the 
estimation of direct and indirect employment in physical terms 
was performed using assumptions based on several sources. The 
total level of full time employment equivalent (FTE) in RES-E 
production was based on National Statistics Data on Environ- 
mental Activities [20]. For the biomass power generation units, 
a number of 220 direct employments were estimated based on a 
survey which allowed concluding that, in average, each power 
generation unit was responsible for the direct employment of 20 
persons [13]. For the solid waste disposal (SWD) units 10 direct 
jobs were considered for each unit, according to a similar study 
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Fig. 1. Cost structure allocation for the production of different RES-E types. 
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Fig. 2. Estimated share of direct and indirect jobs in 2008. 


[31]. For the mini-hydro power generation units and since most of 
these units operate automatically and with few maintenance 
requirements, we have considered 1.5 persons employed for each 
unit (see also ADENE [2]). Big hydro plants only employ perma- 
nently 85 people [28]. For wind we have assumed according to 
EWEA [10] that direct jobs only refer to operation and main- 
tenance which account for 11% of 800. PV is responsible for about 
15 direct jobs in Amareleja plant [1] and an average of 10 in other 
two PV plants [28]. For geothermal energy we have considered 
the same factor of 0.74 jobs per MW installed used in [15] and in 
case of ocean energy we have considered a factor of 0.32 jobs per 
MW installed used in [15,9]. The remaining employments were 
allocated to biogas. 

New I-O vectors have been constructed for each source: hydro, 
wind, PV, biomass (SWD, biomass and biogas) and geothermal. 
Fig. 1 shows the intermediary inputs structure used to build these 
new I-O vectors, which was based on [29] and on the Portuguese 
RES-E production share. 

Fig. 2 illustrates the estimated share of direct and indirect jobs 
for each RES-E type on total direct and indirect jobs, respectively. 
Fig. 3 shows the share of direct and indirect jobs for each RES-E 
type on each RES-E type total jobs. 

The upper bound imposed on RES-E for the time horizon of 
2020 was based on the RES-E production capacity by the 
Portuguese National RES-E Action Plan and considering the same 
rate of production of 2010. We have also considered an upper 
target on the level of FTE consistent with the values attained 
in 2010. 


Geothermaland others 


m Direct jobs 


m Indirect jobs 


Hydro 


Biomass 


0% 20% 40% 60% 80% 100% 


Fig. 3. Share of direct and indirect jobs for each RES-E type in 2008. 


Interval coefficients were only considered for the economic 
constraints and environmental defining constraints. Different 
scenarios were also considered for RES-E in the framework of 
these interval coefficients. 


6. Some illustrative results 


The first solutions were obtained considering a business as 
usual scenario (BAU) corresponding to the base year of this study 
(2008) and three distinct scenarios for the share of RES-E regard- 
ing their corresponding shares on total RES-E production (see 
Table 1). Finally, a scenario was also considered where the RES-E 
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basic equipment (either from rubber and plastics, basic metals, 
metal products, electrical equipment or machinery sectors) is 
domestically produced (see Fig. 10). 


Table 1 
Share of each RES-E on total RES-E production. 


BAU (2008) (%) SCENI (%) SCENII (%) SCENIII (%) 


Hydro 39 45 50 25 

Biomass 30 20 20 35 

Wind 30 30 28 30 

PV 0,2 3 1 9 

Geothermal and others 1 2 1 1 
Table 2 


Extended pay-off tables for each RES-E scenario. 


449 


6.1. Different shares of RES-E 


An extended pay-off table containing all individual optimal 
values for each objective function either with the broadest 
versions of the feasible region and most favourable versions of 
the objective functions (best case scenario) or with stringent 
versions of the feasible region and less favourable versions of the 
objective functions (worst case scenario) is presented for each 
RES-E scenario considered (see Table 2). The main diagonal of 
each pay-off table corresponds to the ideal solution in the best 
and worst case scenarios considered (the individual optimal 
solutions for each objective function either in a best case or worst 
case scenario). GDP (gdp=[gdpworst, dPpbest]) is given in millions 
of Euros, the level of Employment (emp=[eMPworst, EMPpest]) is 


BAU Max gdp Max emp Max res-e SCEN I Max gdp Max emp Max res-e 
Sol 1 Sol 2 Sol 3 Sol 4 Sol 5 Sol 6 Sol 7 Sol 8 Sol 9 Sol 10 Sol 11 Sol 12 
gdpPegest 236.144 216,318 222,229 203,415 202,259 177,510 gdPsest 236,118 216,318 207,821 176,106 154,200 118,631 
gdPwors 218,492 200,866 205,437 188,441 187,300 163,331 gdPworst 218,492 201,557 194,365 165,943 143,830 112,809 
eMpsgest 5,435 4,602 5,600 5,600 4,520 4,116 eMpsest 5,375 5,039 5,600 5,600 3,434 2,874 
eMPworst 5,435 4,602 5,600 5,600 4,520 4,116 eMPworst 5,375 5,039 5,600 5,600 3,434 2,874 
TeS-egest 1627,764 1954211 1627,764 1941,765 5010,039 5010,039 res-epess 2273,369 1953,950 2015,127 1794410 5010,039 5010,039 
TeS-eworst 1627,764 1954211 1627,764 1941,765 5010,039 5010,039  res-eworss 2273,369 1953,950 2015,127 1794410 5010,039 5010,039 
SCEN II Max GDP Max EMP Max RES-E SCEN Ill Max GDP Max EMP Max RES-E 
Sol 13 Sol 14 Sol 15 Sol 16 Sol 17 Sol 18 Sol 19 Sol 20 Sol 21 Sol 22 Sol 23 Sol 24 
gdpPegest 236,127 216,318 196,093 180,306 123,313 141,612 gdPsest 236,112 216,318 236,009 186,816 138,486 118,535 
gdPworst 218,492 202,254 183,379 169,797 117,341 131,826 gdPworst 218,492 201,524 218,492 174,630 129,971 112,718 
eMpesest 5,395 4,595 5,600 5,600 2,711 3,327 emMpPsest 5,367 5,344 5,600 5,600 3,120 2,872 
eMpwors 5,395 4,595 5,600 5,600 2,711 3,327 eMPworst 5,367 5,344 5,600 5,600 3,120 2,872 
TeS-egest 2276566 1962,095 1829,883 1827,200 5010,039 5010,039 res-egest 1855,286 4782,185 1855,286 1855,286 5010,039 5010,039 
TeS-Cworst 2276,566 1962,095 1829,883 1827,200 5010,039 5010,039  res-eworss 1855,286 4782,185 1855,286 1855,286 5010,039 5010,039 
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Fig. 4. Contribution of each RES-E type to DE and IE in each solution obtained in Table 2. 
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given in thousand of persons and RES-E (res-e=[res-eworst, res- 
€best]) production is given in toe (tons of oil equivalent). 

The extended pay-off tables obtained (see Table 2) highlight 
the antagonist nexus between economic growth and RES-E 
production and the overall employment attained in the economy 
and RES-E production. Whenever RES-E production is maximized 
(solutions 5, 6, 11, 12, 17, 18, 23 and 24) the level of GDP and the 
overall level of employment attained suffer a negative impact. 
This fact confirms the main criticisms on RES-E production, 
usually citing that Government subsidies for their production 
may drive up costs and cost jobs or may furthermore crowd out 
other business investment. 

Fig. 4 illustrates the contribution of each RES-E both to direct 
(DE) and indirect employment (IE) in the economy. The NES 2020 
targets regarding RES-E employment are achieved with the 
maximization of RES-E production in all the scenarios considered 
(see Fig. 5), but compromising economic growth and the global 
level of employment attained. Scenario III with the highest 
increase of PV on RES-E production also allows achieving NES 


2020 targets regarding RES-E employment in solution 20 (with 
the optimization of GDP). 

According to BAU scenario RES-E production of hydro and 
biomass contribute with the highest employment shares when- 
ever GDP and the level of the overall employment are maximized. 
PV achieves the highest share on RES-E employment if total RES-E 
production is maximized according to all scenarios herein 
considered. 

However, neither RES-E job studies nor their critiques typically 
include avoided environmental costs or other potential benefits 
(e.g., less imported fossil fuel). Therefore, the assessment of 
imported primary and secondary energy (in toe) and CO, emis- 
sions (in Giga grams—Gg) resulting from fossil fuel combustion 
either in best or worst case scenarios (with more efficient and less 
efficient carbon intensity factors) are illustrated for each solution 
in Figs. 6 and 7, respectively. 

The antagonism between economic growth and the environ- 
mental impacts attained becomes clear in Fig. 7, with CO2 
emissions reaching the highest levels in the solutions which 
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Fig. 5. Contribution of RES-E to total DE and IE in each solution obtained in Table 2. 
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Fig. 6. Energy imports in each solution obtained in Table 2. 
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Fig. 7. CO2 emissions from fossil fuel combustion in each solution obtained in Table 2. 


Table 3 
Minimization of the lower bound of the worst possible deviation. 


BAU SCEN I 


SCEN II SCEN III 


Optimistic stance Optimistic stance 


Optimistic stance Optimistic stance 


Best scenario Worst scenario Best scenario 


Worst scenario 


Best scenario Worst scenario Best scenario Worst scenario 


Sol 25 Sol 26 Sol 29 Sol 30 Sol 33 Sol 34 Sol 37 Sol 38 
SdPpest 218,152 216,318 233,578 204,511 218,152 202,279 213,678 204,730 
SdPworst 204,040 201,624 218,492 191,201 204,040 189,111 198,956 191,404 
emPgest 5,600 5,600 5,600 5,600 5,600 5,600 5,600 5,600 
eMPworst 5,600 5,600 5,600 5,600 5,600 5,600 5,600 5,600 
TES-Cpest 5010,039 5010,039 5010,039 5010,039 5010,039 5010,039 5010,039 5010,039 
TeS-Eworst 5010,039 5010,039 5010,039 5010,039 5010,039 5010,039 5010,039 5010,039 
Table 4 
Minimization of the upper bound of the worst possible deviation. 
BAU SCEN I SCEN II SCEN III 


Pessimistic stance Pessimistic stance 


Pessimistic stance Pessimistic stance 


Best scenario Worst scenario Best scenario 


Worst scenario 


Best scenario Worst scenario Best scenario Worst scenario 


Sol 27 Sol 28 Sol 31 Sol 32 Sol 35 Sol 36 Sol 39 Sol 40 
SdPpest 218,505 216,318 218,505 216,318 218,505 216,318 218,505 216,318 
SdPworst 204,368 202,258 203,637 201,741 203,630 202,246 203,580 202,249 
eMPBest 5,091 5,600 5,065 4,362 4,457 4,451 4,903 4,682 
eMPworst 5,091 5,600 5,065 4,362 4,457 4,451 4,903 4,682 
TES-Cpest 4599,934 3027,470 3246,139 3591,141 4722,173 3027,470 4828,455 3780,650 
TeS-Cworst 4599,934 3027,470 3246,139 3591,141 4722,173 3027,470 4828,455 3780,650 


optimize GDP either in a best or worst case scenario or with 
different shares of RES-E production. As it would be expected with 
the maximization of RES-E production the lowest levels of CO2 
emissions are always obtained. Secondary energy imports also 
reach the lowest levels with the maximization of RES-E produc- 
tion (see Fig. 6). Nevertheless, the weight of RES-E production in 
gross domestic electricity consumption (production plus imports 
less exports) indicates that the maximization of RES-E production 
will lead to higher levels of electricity production than the 
domestic needs. 

The approach used to obtain compromise solutions to the 
MOLP model based on I-O analysis with interval coefficients 
considers two surrogate deterministic problems, considering the 
minimization of the worst possible deviation of the interval 


objective functions from their corresponding interval ideal solu- 
tions [22]. The interval ideal solutions were computed (see 
Table 2) considering both extreme versions of the objective 
functions and the feasible region. The solutions are then obtained 
by minimizing the upper bound (pessimistic stance) or the lower 
bound (optimistic stance) of the worst possible deviation of each 
objective function from its interval ideal solution. 

The values attained for each objective function in best and 
worst case scenarios for each scenario of RES-E production, 
considering the minimization of the lower bound and the mini- 
mization of the upper bound of the worst possible deviation of 
each objective function from its interval ideal solution, are given 
in Tables 3 and 4, respectively. For illustrative purposes we will 
only describe further the solutions obtained with a pessimistic 
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Fig. 8. Contribution of RES-E to DE and IE in each solution obtained in Table 4. 
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Fig. 9. Energy imports and CO, emissions in each solution obtained in Table 4. 


stance (which in our opinion are closer to prospective reality). 
Solutions 25, 26, 29, 30, 33, 34, 37 and 38 that are obtained with 
an optimistic stance (that is, with the minimization of the lower 
bound of the worst possible deviation) always achieve in every 
RES-E scenario herein considered the optimal level of RES-E 
production and the optimal level of employment (see Table 3). 

The solutions obtained by adopting a pessimistic stance allow 
reaching a weight of RES-E production in gross domestic elec- 
tricity consumption from 57% (solution 31) to 115% (solution 39). 
Solution 28 (Scenario BAU) allows obtaining the highest overall 
employment level; however, the level of total RES-E employment 
(22,531 jobs) obtained in this solution is quite far from its NES 
2020 target (see Fig. 8). Solutions 39 and 40 allow reaching NES 
2020 target, but at the expense of a lower overall employment 
value, regarding solution 28. Nevertheless, solution 39 has the 
third highest overall employment level, regarding these last 
solutions obtained (Table 4). 

Although solution 28 allows achieving the highest overall 
employment, it also leads to the highest levels of primary energy 
imports (Fig. 9) and the second highest levels of CO2 emissions. 
On the other hand, solution 39 (with the highest level of RES-E 
production regarding solutions of Table 4) allows reaching the 
lowest levels of secondary energy imports. 


6.2. RES-E basic equipment is domestically produced. 


A scenario was also considered where the basic RES-E equip- 
ment (either from rubber and plastics, basic metals, metal 
products, electrical equipment or machinery sectors) is consid- 
ered to be domestically produced, thus internalizing the import 
values on domestic production coefficients (see Fig. 10). A pay-off 
table containing all individual optimal values for each objective 
function either in a best or worst case scenario is also presented 
for this additional scenario (see Table 5). 

The pay-off table obtained (see Table 5) highlights once more 
the antagonist nexus between economic growth and RES-E 
production and the overall employment attained in the economy 
and RES-E production (the lowest levels of RES-E production are 
obtained with the maximization of GDP—solution 41) and the 
lowest level of overall employment is obtained in the solution 
which optimizes RES-E production—solution 45). Nevertheless, 
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Fig. 10. Coefficients of production of RES-E (Euro/Euro). 


Table 5 
Pay-off table in the best and worst case scenario. 


SCEN IV Max GDP Max EMP Max RES-E 
Sol41  Sol42  Sol43  Sol44 Sol45  Sol46 

ZdPrest 236,144 216,318 206,991 216,318 233,487 211,533 
gdpworse 218,492 202,257 193,588 200,353 218,492 197,028 
eMPpest 5,436 5,600 5,600 5,600 4,860 5,600 
emPworst 5,436 5,600 5,600 5,600 4,860 5,600 
TeS-Cpest 1627,764 1893,289 1775,402 2031,702 5010,039 5010,039 
res-eworst 1627,764 1893,289 1775,402 2031,702 5010,039 5010,039 


with this new data set the employment level reaches the upper 
bound imposed in solutions 42 (maximum Of GDP in a worst case 
scenario), 43 and 44 (maximum Of Employment) and 46 (maximum Of 
RES-E production in a worst case scenario). With this new 
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scenario of coefficients (considering domestic production of RES-E 
basic equipment) it is possible to reconcile the antagonism nexus 
of RES-E production and economic growth and the overall 
employment attained (see, for example, solution 46). 

Fig. 11 illustrates the contribution of RES-E both to direct 
(DE) and indirect employment (IE) in the economy. The NES 2020 
targets regarding RES-E employment are achieved only with the 
maximization of RES-E production, but without compromising 
significantly economic growth and the overall level of employ- 
ment (see solution 46). 

The assessment of imported primary and secondary energy 
(in toe) and CO, emissions (in Gg) resulting from fossil fuel 
combustion either in best or worst case scenarios (with more 
efficient and less efficient carbon intensity factors) are illustrated 
for each solution in Fig. 12. 

The antagonism between economic growth and the environ- 
mental impacts attained becomes once more clear in Fig. 12, with 
CO- emissions reaching the highest levels in the solutions which 
optimize GDP in a best case scenario (solution 41). With the 
maximization of RES-E production low levels of CO, emissions are 
always obtained, albeit GDP is not much compromised in solu- 
tions 45 and 46. Secondary energy imports also reach the lowest 
levels with the maximization of RES-E production (see solutions 
45 and 46 of Fig. 12). Nevertheless, the weight of RES-E produc- 
tion in gross domestic electricity consumption (production plus 
imports less exports) indicates that the maximization of RES-E 
production will lead once more to higher levels of electricity 
production than the domestic needs. 

The values reached for each objective function in best and 
worst case scenarios, considering the minimization of the upper 
bound of the worst possible deviation of each objective function 
from its interval ideal solution are given in Table 6. 

Surprisingly, the best employment and RES-E output results 
are obtained with an optimistic stance, but at the expense of a 
reduction of GDP, highlighting once more the trade-off between 
RES-E output and GDP. 

Fig. 13 illustrates the contribution of RES-E both to direct and 
indirect employment in the economy. The NES 2020 targets regard- 
ing RES-E employment are achieved with the optimistic stance 
(solutions 47 and 48), mainly because of the PV sector. 
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The assessment of imported primary and secondary energy 
(in toe) and CO, emissions (in Gg) resulting from fossil fuel 
combustion either in best or worst case scenarios (with more 
efficient and less efficient carbon intensity factors) is illustrated 
for each solution in Fig. 14. 

Although with an optimistic stance the output of the RES-E 
sector reaches a higher value, higher CO, emissions and primary 
imports of energy are obtained with this new scenario of 
coefficients. This is the result of considering the RES-E basic 
equipment domestically produced from activity sectors (rubber 
and plastics, basic metals, metal products, electrical equipment or 
machinery sectors) highly intensive on energy consumption. 

The solution search process for prospective purposes might 
continue as long as the scrutiny of new solutions is needed. Other 
scenarios could also be considered in order to compute new solutions. 
In this context, another limitation regarding this type of modelling 
approach refers to the I-O data present in this study which is from 
2008 (nevertheless, the I-O tables used herein were released in 
December 2011), being somehow out-of-date, because of the recent 
Portuguese effort in boosting green investment, namely RES-E pro- 
duction. We have tried to overcome this particular limitation by using 
interval data and distinct scenarios. Nevertheless, it is impossible to 
tackle the entire uncertainty herein involved. 


Table 6 
Minimization of the upper bound and lower bound of the worst possible deviation. 


RES-E basic equipment domestically produced 


Optimistic Stance Pessimistic stance 


Best scenario Worst scenario Best scenario Worst scenario 


Sol 47 Sol 48 Sol 49 Sol 50 
gdPpest 210,019 205,768 218,505 216,318 
gdPworss 196,423 192,378 204,367 202,254 
eMPpest 5,600 5,600 4,614 4,644 
eMPworst 5,600 5,600 4,614 4,644 
reS-egest  5010,039 5010,039 3246,139 3027,470 
res-eworst 5010,039 5010,039 3246,139 3027,470 
SCEN IV 


Fig. 11. Contribution of RES-E to total DE and IE in each solution obtained in Table 5. 
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Fig. 12. Energy imports and CO, emissions in each solution obtained in Table 5. 
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Fig. 14. Energy imports and CO2 emissions in each solution obtained in Table 6. 


7. Conclusions 


A model approach is herein presented which entails a pro- 
spective assessment of several factors with impact on energy 
sustainable systems, incorporating energy, environmental, social 
and economic concerns. The analysis performed highlights the 
trade-offs among three objective functions: maximization of GDP, 
maximization of the overall employment level and maximization 
of RES-E production, obtaining a comprehensive analysis of the 
RES-E industry employment. The prospective results obtained 
according to current available data suggest that the targets 
imposed regarding employment generation in the RES-E sectors 
in the National Energy Strategy for 2020 are achieved only under 
extreme assumptions (that is, with the maximization of RES-E 
production) either according to the several RES-E share scenarios 
considered or even if we account for the entire domestic produc- 
tion of the RES-E basic equipment. The results herein achieved 
also allow concluding that with the maximization of RES-E 
production, energy imports, mainly of secondary energy, might 
be reduced as well as CO2 emissions, showing the positive 
influence of RES-E on energy dependence and on the environ- 
ment. Nevertheless, when the RES-E basic equipment is domes- 
tically produced, the maximization of RES-E production does not 
allow attaining the lower levels of CO emissions or the lower 
levels of primary energy imports as it would be expected without 
further analysis, because of the high energy intensity of the 
activity sectors involved in the production of the inputs of these 
sectors. Although, at a local level RES-E production has the merit 
of CO- emissions reduction, at a national level the conclusion is 
not straightforward. On the other hand, the impact of RES-E on 
economic growth is not direct, since there is a trade-off high- 
lighted with the solutions herein analyzed. In fact, this study 
brings to light the antagonist nexus between economic growth 
and RES-E production on one hand and the overall employment 
attained in the economy and the RES-E production on the other 
hand, particularly if the RES-E basic equipment is not mainly 
domestically produced. This fact indicates that the production of a 
significant part of the manufacturing activities domestically is an 
ideal prerequisite for maximizing the positive socio-economic 
effects of the development of RES-E in the Portuguese economy, 
albeit of some unrealistic nature. 


Work is currently under way in order to encompass other 
impacts in this kind of models, namely the impacts of manufac- 
turing and installing energy efficiency measures and the induced 
effects of energy savings and costs reduction. A thorough atten- 
tion will also be paid to the choice of policies which have the 
greatest benefit to cost ratio and on how economic shifts/ 
dislocations should be minimized. 
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